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Spatial Correlation Function of the Chandra Selected Active 

Galactic Nuclei 

Y. R. F. Mushotzky^, A. J. Barger3'^ L. L. Cowie^ 

ABSTRACT 

We present the spatial correlation function analysis of non-stellar X-ray point 
sources in the Chandra Large Area Synoptic X-ray Survey of Lockman Hole 
Northwest (CLASXS). Our 9 ACIS-I fields cover a contiguous solid angle of 0.4 



>■ ■ deg^ and reach a depth of 3 x 10"^^ ergcm~^s~^in the 2-8 keV band. We 

^ ■ supplement our analysis with data from the Chandra Deep Field North (CDFN). 

^ ' The addition of this field allows better probe of the correlation function at small 

■ scales. A total of 233 and 252 sources with spectroscopic information are used in 
^ '■ the study of the CLASXS and CDFN fields respectively. 

■ We calculate both redshift-space and projected correlation functions in co- 
Q^' moving coordinates, averaged over the redshift range of 0.1 < 2; < 3.0, for both 
O ■ CLASXS and CDFN fields for a standard cosmology with = 0.73, Qm = 0.27, 
^ • and h = 0.71 {Hq = lOOh km s~^ Mpc^^). The correlation function for the 

■ CLASXS field over scales of 3 Mpc< s < 200 Mpc can be modeled as a power- 
law of the form ^(s) = {s/sq)'^, with 7 = l.Gto.t and sq = S.OtH Mpc. The 

^ ■ redshift-space correlation function for CDFN on scales of 1 Mpc< s < 100 Mpc 

■ is found to have a similar correlation length sq = 8.55t.o'Jfl Mpc, but a shallower 
slope (7 = 1.3 ± 0.1). The real-space correlation functions derived from the pro- 
jected correlation functions, are found to be tq = 8.1^2:2 Mpc, and 7 = 2.1 ± 0.5 
for the CLASXS field, and tq = b.StH Mpc, 7 = l-SStoH for the CDFN field. 
By comparing the real- and redshift-space correlation functions in the combined 
CLASXS and CDFN samples, we are able to estimate the redshift distortion 
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parameter (3 — 0A± 0.2 at an effective redshift z — 0.94. We compare the corre- 
lation functions for hard and soft spectra sources in the CLASXS field and find 
no significant difference between the two groups. We have also found that the 
correlation between X-ray luminosity and clustering amplitude is weak, which, 
however, is fully consistent with the expectation using the simplest relations be- 
tween X-ray luminosity, blackhole mass, and dark halo mass. 

We study the evolution of the AGN clustering by dividing the samples into 
4 redshift bins over 0.1 Mpc< z <3.0 Mpc. We find a very mild evolution in 
the clustering amplitude, which show the same evolution trend found in optically 
selected quasars in the 2dF survey. We estimate the evolution of the bias, and 
find that the bias increases rapidly with redshift {b{z = 0.45) = 0.95 ± 0.15and 
b{z = 2.07) = 3.03 ± 0.83). The typical mass of the dark matter halo derived 
from the bias estimates show little change with redshift. The average halo mass 
is found to be log {Mhaio/Mo) --12.1. 

Subject headings: cosmology: observations — large-scale structure of the universe 
— x-rays: diffuse background — galaxies: nuclei 

1. Introduction 

Structure formation and evolution in the universe and the formation and growth of 
supermassive black holes (SMBHs) are two fundamental problems in astronomy which are 
still not fully understood. While recent progresses in the cosmic microwave background, 
the high redshift type la supernovae survey, and the large optical surveys have significantly 
improved our understanding of the evolution of large scale structure, there are still several 
gaps in the picture of structure formation. The data at redshift of ~ 1, where most of the 
cosmic star formation might have taken place, is still very limited. On scales of galaxies 
and cluster of galaxies, the feed back process from galaxies or AGNs could significantly 
alter structure formation models where gravitation is the only driving force. The clustering 
of active galactic nuclei (AGNs) provides unique path to the solution of these problems 
because (1) the AGNs are often bright compared to normal galaxies and are easily seen 
at large cosmological distance; (2) AGNs trace the violent growth phase of SMBHs and 
hence their clustering properties provide a link between the dark matter halo to the AGN 
activity. Large scale AGN surveys have been traditionally carried out in the optical band 
with dedicated telescopes. The most recent of these are the Sloan Digital Sky Survey (SDSS, 
Schneider et al. 2004) and the Two Degree Field Survey (2dF, Groom et al. 2005, C05 
hereafter). These surveys have demonstrated that the clustering of AGNs can be used to 
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measure cosmological parameters (Croom et al. 2004; Outram et al. 2004), and to constrain 
gravitational lensing (Myers et al. 2003). However, as has been found in recent Chandra 
and XMM-Newton deep surveys, a large fraction of X-ray detected AGNs show little or no 
activity in optical observations (e.g. Barger et al. 2005), most probably due to obscuration 
(Fabian & Iwasawa 1999) but possibly with some contribution from light dilution of the host 
galaxy (Moran, Filippenko, & Chornock 2002), though Barger et al. (2005) argue this is 
not a dominant factor. This results in a large fraction of AGNs being missed in the optical 
surveys. On the other hand, hard X-rays (> 2 keV) is almost unaffected by obscuring column 
densities Nh < 10^^ cm~^, and the X-ray emission from the host galaxies is low compared to 
the AGNs. Thus hard X-rays are at present the best energy band to find AGNs (Mushotzky 
2004). Recent optical follow-ups of Chandra deep surveys have revealed that the hard X- 
ray sources are mostly found around z ~ 1 instead of 2; ~ 2 as seen in optically selected 
quasar samples. The low redshift population is dominated by non-broadline objects, while 
broadline AGNs are found mostly at higher redshifts (StcfTcn et al. 2003). Given these 
new discoveries, it is important to know how the clustering properties of X-ray and optical 
selected AGN differ. 

The most extensive X-ray AGN surveys so far performed used the ROSAT telescope 
(Mullis et al. 2004). Because the telescope is not sensitive above 2 keV, ROSAT misses 
a large fraction of hard X-ray sources . The relatively poor spatial resolution of ROSAT 
also hmits the accuracy of optical identifications. Most of the ROSAT detected AGNs 
show broad emission lines in their optical spectra. Both the optical quasar surveys and the 
ROSAT surveys suggest that AGNs have correlation properties similar to the local galaxies. 
The results seem to be indcpcndant of the sample medium redshifts. This result is puzzling 
because AGNs are believed to be preferentially form in high density peaks where interactions 
between galaxies are more common, and interactions in turn arc thought to be crucial in 
AGN fueling (Di Mattco, Springel, & Hernquist 2005). The mass of the dark matter halos 
that host AGNs are hence likely to be more massive. 

The clustering results on hard X-ray AGNs are so far contradictory. Earlier studies of 
a small number of individual Chandra fields seem to indicate that the hard band number 
counts in these small fields has fiuctuations larger than expected from Poisson noise (Gowie 
et al. 2002; Manners et al. 2003) but the result is contradicted with larger samples of 
Chandra fields (Kim et al. 2004). Basilakos et al. (2004) found a 4(7 clustering signal in 
hard X-ray sources at f2-8keV > 10~^^ erg cm~^ s"^ using angular correlation functions on a 
XMM detected AGN sample from a 2 deg^ survey. A similar result was also found earlier in 
our 0.4 deg^ Chandra Held (see below) using the count-in-cells technique (Yang et al. 2003). 
Using the Limber equation Basilacos et al. (2004) argue that the hard X-ray sources are 
likely to be more strongly clustered than the optically selected AGNs. Gilli et al. (2003) 
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reported the detection of large angular-redshift clustering in the Chandra Deep field South, 
which seems to be dominated by hard X-ray sources. Using the projected correlation function 
for the optically identified X-ray sources from the Chandra Deep field North (CDFN) and 
South (CDFS), Gilli et al. (2005) found that the average correlation amplitude in the CDFS 
is higher than that in the CDFN, and the latter is consistent with the correlation amplitude 
found in optically detected quasars. 

In this paper, we report our spatial correlation function analysis of the optically iden- 
tified X-ray sources in the Chandra Large Area Synoptic X-ray Survey of the ISO Lockman 
Hole Northwest region (CLASXS). CLASXS is so far the largest contiguous Chandra deep 
field with a high level of spectroscopic identifications. The size of the field is chosen to 
reduce the cosmic variance in the X-ray background to ~ 10% (Yang et al. 2004). For 
comparison, we have also analyzed the correlation functions for the CDFN field, using the 
pubhshed X-ray catalog by Alexander et al. (2003) and the optical catalog of Barger et al. 
(2003). Because the two surveys use basically the same optical instruments in the follow-up 
observations, and thus have the same accuracy in redshift measurements, the comparison is 
relatively straight forward. The LogN-LogS of the CDFN agrees well with that of CLASXS, 
also indicating that the CDFN is a "typical" field. The depth of the CDFN is very useful in 
probing the correlation function at small separations. In § 2 we summarize our observations 
and data analysis. In § 3 we discuss the methodology we use in the clustering analysis. The 
results of the correlation functions are presented in § 4. The evolution of AGN clustering is 
presented in § 5. In § 6 we discussion the implications of our results. Finally summarize our 
results in § 7. Throughout this paper, unless noted otherwise, we assume Hq — 71 and a flat 
universe with = 0.27 and = 0.73. 



2. Observations and data 

CLASXS is a 0.4 deg^ contiguous held centered at a = 10'*34'", S = 57°40' (J2000) in 
the very low galactic absorption Lockman Hole Northwest region. It is the deepest 170//m 
survey held observed by ISOPHOT instrument on board ISO, and has recently been observed 
by the Spitzer Space telescope (Lonsdale et al. 2004). The Chandra observation consists of 
9 ACIS-I pointings separated from each other by ~ 10' to allow a close to uniform sky 
coverage. The center field has an exposure time of ~ 70 ks while the other eight pointings, 
have exposure times of ~ 40 ks. The exposure time were designed to give a uniform fiux 
limit of /2-8 kev ~ 5 X 10~^^ erg cm~^ s^^ , which is about a factor of 2 below the "knee" 
of the LogN-LogS curve. This choice of sensitivity allows a proper sampling of the X-ray 
background sources and also achieves a highest source finding efficiency. 
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The sub-arcsecond spatial resolution of Chandra observatory allows an unambiguous 
optical identification of the X-ray sources, particularly, for those which appear to be normal 
galaxies in optical band. Combined with follow-ups using the large Keck and Subaru optical 
telescopes, we identified and measured the redshifts in a large fraction of the X-ray detected 
AGNs in our survey. The details and the catalogs of the survey can be found in Yang et al. 
(2004) and Steffen et al. (2004). We performed spectroscopic observations for ~ 90% of 
the 525 detected X-ray sources. A total of 272 spectroscopic redshifts have been obtained, 
while the spectra of the rest of the sources have a signal-to-noise ratio too low to obtain 
secure redshift measurements. The redshift distribution of the identified sources are shown 
in Figure 1. The fraction of sources with spectroscopic redshift as a function of hard X-ray 
fiux is shown in Figure 2. 

The 2 Ms CDFN is so far the deepest Chandra field, reaching a fiux limit of f2~skev ~ 
1.4 X 10"^^ erg cm~^s~^ (Alexander et al. 2003). This is ~ 20 times deeper than the 
CLASXS field. The areal density of sources in CDFN is also ~ 5 times higher. The optical 
observation were performed using the same telescope as CLASXS (Barger et al. 2003). 
We use the published catalog, which contains 306 sources with spectroscopic redshift. The 
redshift distribution of the CDFN sources is also shown in Figure 1. The fainter X-ray 
sources in the CDFN are more likely to be found at low redshift, z < 1, compared to the 
CLASXS sources. 



To quantify spatial clustering in a point process, the most commonly used technique is 
the two point correlation function. In short, a two point correlation function measures the 
excess probability of finding a pair of objects as a function of pair separation (Peebles 1980). 



where no is the mean density and r is the comoving distance between two sources. 

Observations of low redshift galaxies and clusters of galaxies show that the correlation 
function of these objects over a wide range of scales can be described by a power-law 
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Methods 



dP^nl[l + i{r)\dVidV2 



(1) 
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(2) 



with 7 ~ 1.6 — 1.9 (Peebles 1980; Peacock 1999). It should be noted that the correlation 
function is in fact a function of redshift, which we will discuss in § 5. Because of the small 
sample sizes of most of the AGN surveys, correlation functions over very wide redshift ranges 
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are commonly used. This only makes sense if the clustering is almost constant in comoving 
coordinates. Fortunately, this is very close to the truth, as we shall see in § 5. 



3.1. Redshift- and real-space Correlation functions 

The nominal distance between sources calculated using the sky coordinates of the sources 
and their redshifts is sometimes called distance in redshift- space, we shall use s instead of 
r to indicate the distance calculated this way. It is apparent that the line-of-sight pecuhar 
velocity of the sources could also contribute to the measured redshift (redshift distortion). 
This effect is most important at separations smaller than the correlation length. The pro- 
jected correlation function, which computes the integrated correlation function along the 
line-of-sight and is not affected by redshift distortion, is often used to obtain the real-space 
correlation function (Peebles 1980). The projection, however, could make the correlation 
signal more difficult to measure. In small fields like the CDFN, the projected correlation 
function is also restricted by the field size, and could be affected by cosmic variance. We 
will calculate both the redshift-space and projected correlation functions in this paper. This 
allows us to estimate the effects of redshift distortion. 

Following Davis & Peebles 1983, we define Vi and V2 to be the positions of two sources 
in the redshift-space, s = Vi — V2 to be the redshift-space separation, and 1 = (vi -|- V2)/2 to 
be the mean distance to the pair of sources. We can then compute the correlation function 
{(rp,7r) on a two dimensional grid, where tt and Vp are separations along and across the 
line-of-sight respectively: 

Tp = Vs ■ s — vr^. (4) 
The projected correlation function is defined as the line-of-sight integration of ^(rp,7r): 

d'Ki{rp,Ti)=i dy^i^r^ + y^), (5) 

where y is the linc-of-sight separation. It has been shown (Davis & Peebles 1983) that, when 
T^max 00, Wp{rp) satisfies a simple relation with the real-space correlation function. If a 
power-law form in Equation 2 is assumed, then 

-.(r,)=r,^-J • (6) 

In practice, the integration is not performed to very large separations because the major 
contribution to the projected signal comes from separations of a few times the correlation 
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length So- Integrating to larger tt will only add noise to the results. After testing various 
scales, we found T^max = 20 — 40 Mpc produces consistent results for our samples. 



3.2. Correlation function Estimator 

To obtain an unbiased estimate of the correlation function, we must correct for selec- 
tion effects. Usually, these selection effects are treated using random samples generated with 
computer simulations. By comparing the simulated and observed pair distributions, the se- 
lection functions effectively cancel. We compute the correlation function using the minimum 
variance estimator 

^ DD - 2DR + RR 
RR 

where DD, DR and RR are the numbers of data-data, data-random and random-random 
pairs respectively, with comoving distances so — A.s/2 < s < so + As/2 (L-S estimator, Landy 
& Szalay 1993). The random catalog is produced through simulations described below to 
account for the selection effects in observations. The random catalog usually contains a very 
large number of objects so that the Poisson noise introduced is negligible. We have checked 
our results using both L-S and the Davis-Peebles estimators (Davis & Peebles 1983) and 
found very good agreement between the two methods. 



3.3. Uncertainties of correlation functions 

There are two terms in the uncertainty of the correlation function: the statistical fluc- 
tuations and the cosmic variance. The statistical uncertainty of the correlation function is 
estimated assuming the error of the DR and RR pairs are zero, and the uncertainty of DD 
is Poissonian, 

In the case of small DD, where y/DD underestimates the error, we use the approximation 
formula (Gehrels 1986) to calculate the Poisson upper and lower hmits. Since the DDs are 
in fact correlated, the use of Poisson errors could underestimate the real uncertainty. In the 
literature bootstrap resampling (Efron 1982) is often used to calculate the errors of the corre- 
lation function. The method is particularly useful in cases when the probability distribution 
function (PDF) of the variable is unknown, or in cases when the variables arc derived from 
Poissonian distributed data using complex transformations, which results in rather complex 
PDFs. Mo, Jing, & Boerner (1992) showed that in the case of large DD, the bootstrap error 
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is ~ -\/3 of the Poisson error. We use Poisson errors in our redshift-space correlation function 

estimates. On the other hand, wc use bootstrap methods when estimating the uncertain- 
ties of the projected correlation function. This is because the numerical integration used in 
Equation 5 make it difficult to apply Poisson errors directly. 

Cosmic variance is known to affect the estimation of the mean density when applied 
to small samples of normal galaxies from optical surveys. Such effect, however, is likely to 
be small on our X-ray selected AGN sample for the following reasons. The volumn of our 
survey is very large compared to the typical pensil-beam optical survey of normal galaxies 
that typically covers very narrow rcdshift ranges. On the other hand, the number density 
of AGNs is much lower than that of the normal galaxies, making it hard to trace individual 
structures at high enough sampling rate. The window function of the spectroscopic follow-up 
in our survey is also very flat over a wide redshift range (except in the redshift desert at 
2; ~ 1.2 — 2). The combination of these factors makes it very difficult for a small number of 
structures been over sampled and thus producing incorrect estimation of the mean density. 
However, for ultra deep surveys with field size of a single Chandra field, small number of 
velocity spikes can indeed affect the correlation analysis, as seen in the case of the Chandra 
Deep Field South. Such structure, however, will afi^ect number counts in the field at fiux 
levels comparable to the depth of CLASXS. Based on the very good agreement among the 
number counts found in the CLASXS, CDFN, and other deep surveys (Yang et al. 2004), 
we believe the uncertainty from cosmic variance on the whole sample is likely to be small. 
However, the cosmic variance effect on subsamples could still be important, as seen in §5. 
In such cases, using statistical uncertainty alone could underestimates the true uncertainty. 



3.4. The mock catalog 

To account for the observational selection and edge effects, we perform extensive simu- 
lations to construct a mock catalog. 

The Chandra detection sensitivity is not uniform because of vignetting effects, quantum 
efficiency changes across the field and the broadening of the point spread functions. The 
consequence is that the sensitivity of source detection drops monotonically with off-axis 
angles. To quantify this we generate simulated observations of our 40 ks and 70 ks exposure 
in both soft and hard bands. Using wavdetect (Freeman et al. 2002) on these images we 
obtain an estimate of the detection probability function at different fluxes and off-axis angles 
(Figure 3). 

With this probability, we can generate randomly distributed sources with the X-ray 
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selection effects to tfie first order. We use tfiis metfiod instead of running detections on a 
large number of simulated images because the detection program runs very slowly on these 
images. We generate source fluxes based on the best fit LogN-LogS from Yang et al. (2004) 
and then "detections" are run on each of the images. The resulting catalogs from all the 
nine simulated images are then merged in the same way as for the real data. The resulting 
random source distribution and the resulting cumulative counts are shown in Figure 4. 

We next consider the optical selection effects. Since our spectroscopic observation is 
close to complete for all sources with R< 24.5, the sky coverage is uniform and only a very 
small number of sources which are very close to each other could be missed. This might 
reduce the power at very small scales. The redsliift distribution of the sources shows a 
very weak dependence on the X-ray flux (Figure 5), which is due largely to the very broad 
luminosity function of AGNs. We can thus "scramble" the observed redshifts and assign 
them to the simulated sample without introducing a significant bias. The major selection 
effect in our optical observation is that the optical identifications are biased toward brighter 
sources. We select the simulated sources based on their X-ray fiux using the best-fit curve 
in Figure 2 as a probability function. The optical selection removes a large fraction of X- 
ray dim sources and therefore reduces the non-uniformity in the angular distribution caused 
by the X-ray selection effects. The redshift of the random sources were sampled from a 
Gaussian smoothed {a^ — 0.2) redshift distribution from the observations. The purpose of 
the smoothing is to remove possible redshift clustering in the random sample but still preserve 
the effect of the selection function. We tested different smoothing scales Az = 0.1 — 0.3 and 
found the resulting correlation function effectively unchanged. 



4. Results 

4.1. Redshift-space correlation function 

We calculate the redshift-space correlation function in the CLASXS sample for non- 
stellar sources with 0.1 < z < 3 and 2-8 keV fiuxes > 5 x 10^^^ergcm~^s~^ , assuming 
constant clustering in comoving coordinates. The total number of sources in the sample is 
233. The median redshift of the sample is 1.2. By comparing the number of detected pairs 
with separations < 20 Mpc with that expected by simulation, we found that on scales of 
20 Mpc, the significance of clustering is 6.7a. 

We use the maximum likelihood method in searching for the best-fit parameters (Cash 
1979; Popowski ct al. 1998; Mullis et al. 2004). The method is preferable to the commonly 
used method because it is less affected by arbitrary binning. The method uses very small 



bins so that each bin contains only 1 or DD pair. In this hmit, the probabihty associated 
with each bin is close to independent. The expected number of DD pairs in each bin is 
calculated using the DR, RR pairs using the mock catalog. The likelihood is defined as 
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where /Xj is the expected number of pairs in each bin and Xi is the observed number of pairs. 
The likelihood ratio defined as 



where Cq is the maximum likelihood. The resulting S approaches the usual distribution. 
Since the maximum-likelihood does not provide a measure of the "goodncss-of-fit", we quote 
the derived from the binned correlation function (as shown in the figures) and the best-fit 
parameters from maximum-likelihood estimates. 

We fit the correlation functions over three separation ranges. In Figure 6 we show the 
correlation function and the best-fit with 3 Mpc< s <200 Mpc. The best-fit parameters for 
all three separation ranges are hsted in Table 1. It is noticeable that the rather large 
seems to suggest that the single power-law model may not be a proper description of the 
data. 

For comparison, we also computed the correlation function of the X-ray sources in CDFN 
in the same redshift interval. We created a mock catalog 50 times larger than the observation. 
The positions and redshifts of the random sources are generated by randomizing the observed 
positions and redshifts. A large Poisson noise was added to avoid artificial clustering in the 
mock catalog. Such randomization is justified because the clustering signal in a small field like 
the CDFN mainly comes from clustering along the line-of-sight direction. The randomized 
sky coordinates are filtered using an image mask to take into account the edge effects. We 
include all the non-stellar sources in the same redshift interval as we use for CLASXS, which 
results in 252 sources in the sample. The best-fit parameters for CDFN field over three scale 
ranges are also shown in Table 1. The correlation function over 1 Mpc< s <100 Mpc is 
shown in Figure 7. 

There is a good agreement of the correlation lengths obtained in the two fields. There 
seems to be a systematic fiattening of the slope at small separations (s ~ 10 Mpc) in 
both samples. When the correlation functions are fitted at small and large separations 
independently, the resulting x^s are systematically smaller. 



S=-2\n{C/Co) 



(10) 
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4.2. Projected correlation function 

The projected correlation function is computed using the methods described in § 3.1. 
To test the method, we first compute the projected correlation function for the CDFN and 
compare the results with that published in Gilli et al. (2005). We selected the same redshift 
interval for the CLASXS field. A two dimensional correlation function is calculated on a 5 x 10 
grid on the (rp,7r) plane. The 5 intervals along Vp axis covers 0.16-20 Mpc. We integrate the 
resulting two dimensional correlation function along the line-of-sight to a TTmax — 20 Mpc. 
Our projected correlation function for CDFN is shown Figure 8, and it agrees perfectly with 
that reported in Gilli et al. (2005) for 2; = — 4, validating the techniques used in this paper. 

We next compute the projected correlation function for the CLASXS field. The corre- 
lation function is calculated on scales of Tj, = 1 — 30 Mpc. The 2-D correlation function is 

integrated to iTmax = 30 Mpc. The result is also shown in Figure 8. The correlation functions 
of the CDFN and CLASXS fields agree in general at Vp ~ 10 Mpc, where both surveys have 
very good S/N. The slope, however, appears to be fiatter in the CDFN field. 

We perform a fit to the correlation functions using Equation 6. The best-fit parame- 
ters for the CDFN are tq = 5.8t\-^^ Mpc, 7 = 1. 38^0:14, and the reduced x^/dof = 2.5/3. This 
is in good agreement with the result from Gilli et al. (2005, Tq = 5.7 Mpc, 7 = 1.42). The 
quoted errors in that paper is smaller than we obtained, but since we adopt a bootstrap error 
instead of Poisson error in this analysis, the difference is expected. The best-fit parameters 
for the CLASXS field are ro = S.ltH Mpc, 7 = 2.ll°J, and the reduced x^/dof = 1.6/4. 
The correlation length appears to be higher than that of the CDFN, but agrees within the 
errors. The slope also seems steeper than that of the CDFN and agrees better with the slope 
of the redshift-space correlation function at > 10 Mpc. Since the CLASXS sample does 
not cover separations < 10 Mpc very well, it is hard to see a slope change in this sample 
alone. Since the CDFN and CLASXS connect very well at separations where both surveys 
are sensitive, we try to model the combined data points with a single power-law. This yields 
r-Q = 6.lt?'o Mpc, 7o = 1.47tg5oj X^/dof = 10.7/9. The reduced x"^ is much worse than 
the two samples fitted separately, but does not reject the hypothesis of a single power- law fit. 
This seems to suggest that the slope of the correlation function flattens at small separations. 
Such a trend need to be looked more closely with future large AGN surveys with better 
signal-to-noise. 
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4.3. 



Redshift distortion 



Rcdshift distortion affects the correlation function (power-spectrum) by increasing tlie 
redsfiift-space correlation amplitude and changing the shape of the 2-D redshift-space corre- 
lation function at small scales (such as the well known "finger-of-God" effect, e.g. Hamilton 
1992). Since our data is too noisy at small separations to detect the effect, we only discuss 
the effect of the amplitude boosting of correlation function in redshift-space. Kaiser (1987) 
showed that to the first order. 



where j3 ~ VIm{z)^'^ /i>{z) and h{z) is bias. In principle, the redshift-space distortion can 
be estimated by comparing ^{s) and i{r). To quantify the effect, we use the correlation 
function estimates at scales where both projected and redshift-space correlation functions 
are well determined. For the CDFN, we chose the correlation function estimates at 10 Mpc 
and find ^(s = 10 Mpc)/^{r = 10 Mpc) = 1.75 ± 0.55, if the best-fit of ^(s) on 1-100 Mpc 
is used. The choice of this scale is justified given that the slope possibly changes below 
and above 10 Mpc, as seen in the projected correlation function. Since the slope of the 
redshift- and real-space correlation function is very similar in the CDFN, the ratio is almost 
constant. For the CLASXS field, we chose to estimate the ratio at 20 Mpc, where the S/N 
is the best. We find ^{s = 20 Mpc)/i{r = 20 Mpc) = 1.73 ± 0.42 by using the best-fit 
on 1-100 Mpc for ^{s). The ratio changes slowly with the scales probed, but is within the 
errors. We find a general agreement between CLASXS and CDFN. To avoid the arbitrary 
choice of scales, and to make the best use of the data, we combine the two samples to study 
the redshift distortion effect on ^(rp,7r). Since the projected correlation function of CDFN 
and and CLASXS agree in general, we are encouraged to assume that the the two samples, 
even with the vast difference in fiux limits, generally trace the large scale structure in the 
same way. 

In Figure 9. we show the combined ^(rp,7r). The contours show no significant signature 
of nonlinear redshift distortion, such as the "finger-of-god" . We fit ^{vp, n) with Equation 11, 
assuming the best-fit parameters for the real-space correlation function from the combined 
sample {tq = 6.1 Mpc, 70 = 1.47), and ignoring the higher order redshift distortions. We 
generate the 2-D correlation function at each grid point. By minimizing by changing /?, 
we found the best-fit (3 — 0.4 ± 0.2, corresponding to C{s)/^{r) ~ 1.3, which agrees with the 
estimates from individual fields above. By fixing = 0.27, we can estimate the bias factor 
of X-ray selected AGNs from (3. The median redshift of the combined sample is 0.94, and 



nM{z = 0) = 0.27 gives Qniz = 0.94) = 0.73. This yields b ^ 2.04 ± 1.02 using the relation 



(11) 



(3 ^ n%yb. 
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4.4. X-ray color dependence 

We further test if there is any differences in clustering properties between the hard 
and soft spectra sources in the CLASXS sample. We use the hardness ratio, defined as 
HR = C2-8keV / Co,5-2keV (where C is the count rate), to quantify the spectral shape of the 
X-ray sources. Correlation functions of soft {HR < 0.7) and hard {HR > 0.7) sources are 
calculated the same way as above. The fraction of broad-hne AGNs is 56.4% in the soft 
sample and 15.4% in hard sample. The median redshifts are 1.25 and 0.94 for soft and 
hard samples, respectively. We compute ^{s) for both soft and hard sources over scales of 
3-200 Mpc. 

Using a maximum- likelihood fit, we found So = 9.61^:^ Mpc, 7 = IMtoi for hard sources 
and So = S-S^lo Mpc, 7 = l-6lo;5 fo^ sources. We found no significant difference in clus- 
tering between the soft and hard sources. This agrees with the results of Gilli et al. (2005). 
It is noticeable that the soft sources have a higher median redshift than the hard sources. 
The interpretation of this result must include evolution effects. To avoid this complication, 
we restricted the redshift range to z = 0.1 — 1.5. The best-fit parameters are Sq = 9.5l|^ Mpc 
{6.2^.41 Mpc) and 7 = l.Tlgie C^-^-o^) fo'^ hard (soft) sources. The difference in clustering 
parameters between soft and hard sources are well within the measurement error. The same 
analysis on CDFN yields similar results. Thus there is no significant dependence of clustering 
on the X-ray color. 

4.5. Luminosity dependence 

The cold dark matter (CDM) model of hierarchical structure formation predicts that 
massive (and hence luminous) galaxies are formed in rare peaks, and therefore should be 
more strongly clustered. This is seen in normal galaxies (e.g. Giavalisco & Dickinson 2001). 
Whether this relation can be extended to X-ray luminosity of AGNS is unknown. This is 
because the X-ray luminosity relates to the dark matter halo mass in a more complex and not 
well understood way. The X-ray luminosity is directly linked to the accretion process, and 
the process is affected by factors such as accretion rate, radiative efficiency, blackhole mass 
and the details of the accretion process. We have shown that at least in broadline AGNs, 
where the blackhole mass can be inferred from the line-width and nuclear luminosity, the 
Eddington ratio is close to constant over two decades of 2-8 keV luminosity (Barger et al. 
2005). If this is the case for all X-ray selected AGNs, we should expect the AGN luminosity to 
be mainly determined by the blackhole mass, which in turn, should be closely related to the 
halo mass (Ferrarese 2002), even though the exact form of this relation is highly uncertain. 
However, optical quasar surveys such as 2dF found little evidence of a correlation between 
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clustering amplitude and ensemble luminosity (COS), perhaps due to the small dynamical 
range in luminosity these surveys probe. The X-ray luminosity of sources in the CLASXS 
and CDFN cover a luminosity range of four orders of magnitudes, making it possible to make 
such a test. 

The 2-8 keV rest frame luminosity is calculated from the hard band fluxes, with a 
K-correction made assuming a power-law spectra with photon index F = 1.8. This yields 



where Lq is the luminosity in observer's 2-8 keV band. In Figure 10 we show vs. redshift 
for both CLASXS and CDFN. For a better comparison of the correlation amplitude, we 
adopt the averaged correlation function within 20/i~^ Mpc, 



The quantity is chosen rather than sq because it measures the clustering (directly linked to 
the rms fluctuations) regardless of the shape of the correlation function. On scales of 20 Mpc 
the clustering is in the hnear regime of density fluctuations. The error in ^(20) is from the 
single parameter la confldence interval obtained by flxing the slope of the correlation function 
to the best-flt. 

We split the CLASXS (CDFN) sample into two subsamples at = 4.5 x 10^^ erg s"^ (3.2 x 
10^^ ergs~^) so that each subsample contain similar number of objects. In Table 2 we 
show the maximum-hkelihood flts as well as C(20)s. It should be noted that the correla- 
tion amplitude is biased in redshift space. The dominant part of this bias is characterized 
in Equation 11. Comparing with other observations (e.g. da Angela et al. 2005), (3 is 
likely a weak function of redshift in the redshift range probed by our sample, with (3 ~ 0.4 
(§4.3), this translates to ^{s)/^{r) ~ 1.3. We correct the ,^(20) 's for this bias by dividing 
them by 1.3. The correlation amplitude for the more luminous sources appears to be higher 
than that of the less luminous sources, which qualitatively agrees with expectations that 
X-ray luminosity reflects the dark matter halo mass. The correlation amplitude for the more 
luminous subsamples are 2.3(7 and 5.7a" higher than that of the less bright subsample in 
the CLASXS and CDFN fields, respectively. However, since the more luminous subsamples 
also are preferentially found at higher redshifts, the evolution in ^(s) should be taken into 
account. 

To reduce this complication, we restrict ourselves to sources within the redshift range 
of 0.3-1.5, where the evolution effect is relatively small (see also § 5). In Figure 11 we show 
Lx vs. ^(20) for both CLASXS and CDFN. By reducing the redshift range, the difference in 
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correlation amplitude between the brighter and dimmer subsample is significantly reduced, 
in the CDFN sample, to merely la. For the CLASXS field, on the other hand, the correlation 
amplitude for both subsamples do not show significant change. For comparison, we also plot 
in Figure 11 the correlation amplitude from the 2dF survey (C05). The X-ray luminosities 
for the QSOs in the 2dF arc obtained by dividing the bolometric luminosities by 35 (Elvis 
et al. 1994). We perform Spearman's p test for correlations between log and ^. We found 
the correlation coefficient p — 0.8 for X-ray samples, or a corresponding null probability of 
20%, indicating a weak correlation between the two quantities. If the 2dF samples are added, 
we found p — 0.1, and a null probabihty of 17%. This means that with the X-ray sample, 
we have detected a weak correlation between clustering and luminosity. We will discuss this 
in § 6.2. 



5. Evolution of clustering 

Measuring the correlation function over a wide redshift range only makes sense if the 
correlation function is a weak function of redshift. The best measurements of clustering of 
2dF quasars at high redshift show that the correlation function indeed exhibits only mild 
evolution (C05). In this section, we test the evolution of clustering of X-ray selected AGNs 
and compare them with other survey results. 



5.1. Chandra Sample 

We study the evolution of clustering in both CLASXS and CDFN samples, using the 
redshift-space correlation function. The sources are grouped in 4 redshift intervals from 0.1 
to 3. The sizes of the intervals are chosen so that the number of objects in each interval 
is similar in the CLASXS sample. This result in a very wide redshift bin above z = 1.5. 
The correlation functions for the CLASXS, CDFN and CLASXS+CDFN fields are shown in 
Figures 12, 13, and 14, respectively. We group the pair separations in 10 bins in these figures 
to show the shape of the correlation function. In some bins there could be no DD pairs, and 
the correlation function is set to -1 without errors. We model the correlation functions using 
single power-laws and fit the data using the maximum-likelihood method. As we mentioned 
earher, the method is not affected by binning. We found on 3-50 Mpc scales that a single 
power-law provides a good fit to the data except, for the the z — 1.5 — 3 interval in the 
CDFN, where the sample is too sparse and have very few close separation pairs, we use a 
separation range of 5-200 Mpc to obtain the fit. The goodncss-of-fit is quantified with x^- 
In the case where empty bins exist, we increase the bin sizes until no bins are empty before 
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we compute the x^- The results are summarized in Table 3 and the ^(20)s as a function 
of redshift are shown in Figure 15. We have tested fitting the correlation functions over 
different scale ranges, and found no significant difference in the resulting ^(20)s. 

There is only mild evolution seen in both the CLASXS and CDFN fields, in agreement 
with the assumption that clustering is close to constant in comoving coordinates. There are 
some small discrepancies between the CLASXS and the CDFN clustering strength. These 
discrepancies give the sense of the field-to-field uncertainty and could be resolved with future 
larger surveys. At the highest redshift, both samples show an increase trend of clustering, 
but only at the < 2a level. 

5.2. Comparing with other observations 

In Figure 16 we plot ^(20) as a function of redshift for CLASXS, the combined CLASXS 
and CDFN, as well as results from the 2dF (C05), the ROSAT North Galactic Pole Survey 
(NGP, Mulhs et al. 2004), and the Asiago-ESO/RASS QSO survey (AERQS, Grazian et al. 
2004). We did not correct for redshift distortion for observations which uses redshift-space 
correlation function. This leads to overestimates of the real-space correlation amplitude. 

Our correlation functions show clear agreement with the evolution trend found in COS. 
However, as seen in § 4.5, our measured correlation amplitude on average appears similar as 
that of 2dF, even though the average luminosity of latter is much higher than that of the 
X-ray samples. We compare the X-ray luminosities of the CLASXS and CLASXS-I-CDFN 
samples with those of the 2dF in Figure 17. The X-ray luminosities of 2dF quasars are 
obtained the same way as in § 4.5. The luminosity difference between the 2dF sample and 
X-ray samples is the largest at low redshift and decreases at higher redshift. At z > 2, the X- 
ray sample and the 2dF samples have similar median luminosity. The similarity in clustering 
amplitude can be understood in the light of the weak correlation between AGN luminosity 
and clustering amplitude found in § 4.5. However, as we will see in § 6.2, a correlation 
between dark halo mass and the X-ray luminosity predicts a rapid increase of correlation 
function above Lx ~ 10 ergs . Therefore, we should expect the optical sample being 
more clustered than X-ray samples at medium redshifts. Such a trend is not clearly seen in 
our samples. 
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6. Discussion 

6.1. Evolution of Bias and the typical dark matter halo mass 

The bias evolution of optical quasar is extensively discussed in COS. They found that 
the bias increases rapidly with redshift(6 ~ (1 + z^). We will follow these arguments to 
estimate the bias evolution of the X-ray samples. 

On scales of 20 Mpc, the clustering of dark matter and AGNs are both in the linear 
regime, i.e., ^(20) < 1. This allows us to measure the bias, defined as the ratio of rms 
fluctuation of the luminous matter (AGNs in our case) and that of the underlying total 
mass, as a function of redshift by comparing the observed correlation function with the 
linear growth rate. In terms of correlation function, the bias can be written as 

b"^ = flight /^mass- (14) 

The averaged correlation function of mass can be obtained using 

where ^2(7) = 72/[(3 — 7)(4 — 7)(6 — 7)2'''], ag — 0.84 is the rms fluctuation of mass at 
z — Q obtained by WMAP observation (Spergel et al. 2003), and we choose the best-flt 
7 ~ 1.5. D{z) is the linear growth factor, for which we use the approximation formula 
from Carroll, Press, & Turner (1992). The redshift-space distortion is taken into account 
to the first order through Equation 11 and the bias factor is solved numerically. The result 
is shown in Table 4. The estimate of b{z = 1) ~ 2.2 in the combined sample agrees with 
the result from the redshift-space distortion analysis in § 4.3. In Figure 18(a) we show the 
bias estimates for the CDFN and CLASXS+CDFN samples. The best-fit model from COS 
qualitatively agrees with the X-ray results. 

The simplest model for bias evolution is that the AGNs are formed at high redshift, and 
evolve according to the continuity equation (Nusser & Davis 1994; Fry 1996). The model is 
sometimes called the conserving model or the test particle model. By normalizing the bias to 
z — 0, the model can be written as 

b{z)^l + {bo-l)/D{z). (16) 

This model is shown in Figure 18(a) as dash-dotted line. The model produces a bias evolution 
which is shghtly too shallow at high redshifts. The correlation function evolution based on 
this model is also shown in Figure 16, where it underpredicts the observed ^. The model 
predicts a decrease of correlation function at high redshift, which is not true based on our 
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results and that of the 2dF. This imphes that the AGNs observed in the local universe are 
unhkely to have formed at z » 2. 

One of the direct predictions of CDM structure formation scenario is that the bias is 
determined by the dark halo mass. Mo & White (1996) found a simple relation between 
the minimum mass of the dark matter halo and the bias b. By adopting the more general 
formalism by Sheth, Mo, & Tormen (2001) we can compute the "typical" dark halo mass 
of the sample. It should be noted that the method assumes that halos are formed through 
violent collapse or mergers of smaller halos and hence is best applied at large separations, 
where the halo-halo term dominates the correlation function. This requirement is apparently 
satisfied by AGNs. Following Sheth, Mo, & Tormen (2001), 

where u = Sc{z)/a{M, z), a = 0.707, c = 0.6. Sc is the critical overdensity. a{M,z) is the 
rms density fluctuation in the linear density field and evolves as 

a{M,z) = ao{M)D{z), (18) 

where (7q{M) can be obtained from the power spectrum of density perturbation P{k) con- 
volved with a top-hat window function W{k), 

aoiM) = ^ y dkeP{k)\Wik)\^ (19) 

At the scale of interest (~ 10 Mpc), the power spectrum can be approximated with a power- 
law, P{k) oc A;", with — 2 < n < — 1 for CDM type spectrum. Integrating Equation 19 
gives 

^o(M) = ^8(^)-("+^)/^ (20) 

where Mg is the mean mass within 8 Mpc. 

We can then solve Equation 17 for halo mass. The resulting mass is shown in Table 4 
and Figure 18(b). Consistent with what's been found in COS for the 2dF, the halo mass does 
not show any evolution trend with redshift. We found < log{Mhaio/ Mq) >~ 12.11±0.29, 
which is consistent with the 2dF estimates (COS, Grazian et al. 2004). 



6.2. Linking X-ray luminosity and clustering of AGNs 

We have shown that over a very wide range of luminosity, the clustering amplitude of 
AGNs change very httle. This allows us to put useful constrains on the correlations among 
X-ray luminosity, blackhole mass Mbh, and the dark matter halo Mhaio- 
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Using the equivalent width of broad emission hnes as mass estimators, Barger et al. 
(2005) found that the Eddington ratio of broadhne AGNs may be close to constant. Since 
the hard X-ray luminosity is an isotropic indicator of the bolometric luminosity, this implies 
that the blackhole mass is linearly correlated with X-ray luminosity. Barger et al. (2005) 
found that 

where L44 is in units of 10^^ ergs~^ . A similar result is found at low redshift using a 
sample of broadline AGNs with mass estimates based on reverberation mapping relations 
(Kaspi et al. 2000; Yang 2005). The relation, however, is only tested for broadline AGNs. 
Deviations from this relation is also expected at low luminosities since many low luminosity 
AGNs tend to have a low Eddington ratio (Ho 2005). 

Blackhole mass have been shown to correlate with velocity dispersion of the spheroidal 
component of the host galaxies (Gebhardt et al. 2000; Ferrarese & Merritt 2000). This 
leads to a linear correlation between Mbh and the mass of spherical component. This 
relation, however, could be different at high redshift (Akiyama 2005). How these relationships 
translate to the Mbh ~ M^aio relation is also unclear and could likely be nonlinear. Ferrarese 
(2002) showed that Mbh - Mhaio can be modeled with a scaling law 

(22) 



108 Mq 'IO12M0 

with K and A determined by the halo mass profile. 

Combining the above and using Equation 17, we can calculate the correlation amplitude 
as a function of X-ray luminosity. In Figure 11 we show the model expectations compared 
with the observations from CLASXS, CDFN and 2dF. In calculating the bias we have as- 
sumed the nonlinear power-law index n = 3 — 7 in Equation 20 (Peacock 1999), with the 
best fit 7 = 1.5. The three lines represent three different halo profiles discussed in Ferrarese 
(2002). We found that the — C,{'20) relation is in fact dominated by the very nonlin- 
ear relation between halo mass and correlation amplitude. The difference between different 
halo profiles is caused mainly by the normalization k, or the fractional mass of blackhole 
mass, rather than the power- law index A. One of the important model predictions is that 
the correlation between X-ray luminosity and clustering is weak below ~ 10^^ erg s~^ and 
increases rapidly above that. The lack of rapid change of correlation amplitude indicates the 
halo mass of AGN cannot be significantly higher than the corresponding threshold. Under 
the assumed cosmology and bias model, the — C(20) relation based on the weak lensing 
derived halo mass profile (Seljak 2002) is consistent with the data, while the NFW profile 
(Navarro, Frenk, & White 1997) and the isothermal profile predicts a too steep correlation 
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amplitude curve at high luminosity. However, we cannot rule out these profiles as a rea- 
sonable descriptions of the AGN host halo because of the uncertainty in the shape of the 
correlation function and the fractional blackhole mass in dark halos at the rcdshift of our 
sample. In Figure 11 we also mark the model dark halo mass corresponding to the Seljak 
(2002) mass profile. The average correlation amplitude of the combined optical and X-ray 
sample (dotted-line) corresponds to a halo mass of ~ 2 x 10^^ Mq. While the luminosity 
in our sample ranges over five orders of magnitudes, the range of halo mass may be much 
smaller. The 2dF sample has a high luminosity but has a similar or lower average correlation 
amplitude than that of the X-ray samples. It is possible that the optical selection technique 
tends to select sources with a higher Eddington ratio. The correlation amplitude of the 
CDFN sample at ~ 10^^ ergs~^ , on the other hand, is higher than the model predictions. 
This is expected because many AGNs with such luminosities are LINERs which are probably 
accreting with a low radiative efficiency. 

It is now clear that the weak himinosity dependence of AGN clustering is consistent with 
the simplest model based on the observed — Mbh and Mbh — Mhaio relations, given the 
large error in the correlation functions. A large dynamical range in X-ray luminosity, as well 
as better measurements of correlation function, are needed to better quantify this relation. 
The luminosity range of the 2dF survey is too small and the optical selection method also 
likely is biased to high Eddington ratio sources. By increasing our current CLASXS field 
by a factor of a few will be helpful in better determine the luminosity dependence of AGN 
clustering, and to put tighter constrains on AGN hosts. 



6.3. Blackhole mass and the X-ray luminosity evolution 

We look again at the MBH^Mhaio relation in the light of the mass estimates of the 
dark matter halos from Chandra samples. If the Ferrarese (2002) relation is independent of 
redshift, the nearly constant dark halo mass implies little evolution for the blackhole mass. 
On the other hand, strong luminosity evolution is seen since 2; = 1.2 in hard X-ray selected 
AGNs (Barger et al. 2005). This imphes a systematic decrease of the ensemble Eddington 
ratio with cosmic time. Barger et al. (2005) showed that the characteristic luminosity of 
hard X-ray selected AGNs 

L* = Lo(^)«, (23) 

where log(Lo/ ergs~^ ) — 44.11 and a — 3.2 for z < 1.2. Ueda et al. (2003) found a similar 
result with a slightly shallower slope. If the typical blackhole mass does not change with 
redshift, the observed luminosity evolution indicates the ensemble Eddington ratio increase 
by a factor of ~ 10 from z — to z — 1. It is hard to understand such a change of the 
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typical Eddington ratio with redshift. One possibility is that a large number of Compton 
thick AGNs at z ~ 1 are missed in the Chandra surveys (e.g. Worsley et al. 2005), leading 
to the observed strong luminosity evolution. 

Alternatively, instead of MBH-Mhaio being independent of redshift, the Mbr-Vc could 
be unchanged with cosmic time, as suggested by Shields et al. (2003). This is theoretically 
attractive because the feedback regulated growth of blackholes implies a constant Mbh^c 
relation. This implies that MsH^Mhaio is in fact a function of redshift. Wyithe & Loeb (2003) 
proposed a model(WL model here after) showing that the blackhole mass inferred from the 
halo mass increases with redshift. Groom et al. (2005) show that this could lead to a close to 
constant Eddington ratio in the 2dF sample if the optical luminosity is used to compare with 
the derived Mbh- Since the correlation function is only a weak function of luminosity, as 
we have demonstrated in § 4.5, it is better to estimate the evolution of the Eddington ratio 
using the characteristic mass of the blackholes from the WL model, and the characteristic 
luminosity from Equation 23. In Figure 19, we show the derived ensemble Eddington ratio, 
assuming the dark halo mass to be constant and log (< Mhaio/MQ >) ~ 12.11. (we adopt 
the normalization of the WL model so that it matches the prediction of Mbh ^ Mhaio with 
a NWF type of halo profile. However, the choice of this normalization is not crucial). In 
the figure, we see a factor of ~ 2 — 3 change in the ensemble Eddington ratio from z = 
to z = 1.2. This change, however, is smaller than the typical scatter in both the luminosity 
and halo mass. 



6.4. Compcirison with normal galctxies 

We now compare our clustering results with those for normal galaxies. Using the Sloan 
Digital Sky Survey First Data Release, Wake et al. (2004) found that the clustering of 
narrow-hne AGNs in the redshift range 0.055 < z < 0.2, selected using emission-hne flux 
ratios, have the same correlation amplitude as normal galaxies. Our samples are not a very 
good probe at these redshifts, and the best clustering analysis at a comparable redshift for 
normal galaxies is from DEEP2 (Goil et al. 2004). At effective redshift Zeff ~ 1, they 
found ro = 3.19 ± 0.51 h"^ Mpc, and 7 = 1.68 ± 0.07, which translates to ^(20) ~ 0.1. 
The correlation amplitude from GLASXS at ^ = 0.9 is ^(20) ~ 0.13+o:}o- The clustering of 
AGNs in GLASXS field appear the same as the clustering of normal galaxies in DEEP2. A 
higher correlation is found in the GDFN. The difference shows the large uncertainty of our 
correlation function estimates. At higher redshifts, the best estimate for galaxy clustering 
is from the so called "Lyman break galaxies" , named after the technique by which they are 
found. Adelberger et al. (1998) found, at a typical 2; ~ 3, these galaxies tend to have similar 
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correlation function as galaxies in the local universe, indicating they are highly biased tracers 
of the large scale structure. In the ACDM cosmology, these authors found b = 4.0 ±0.7. This 
is very similar to the bias found in the highest bin of our Chandra fields {b = 3.03 ± 0.83), 
which has a median redshift of ~ 2.0. If we extrapolate the bias of the X-ray sources to 
z = 3, the bias of X-ray sources should be ~ 4 — 5, consistent with the clustering strength 
of Lymann break galaxies. 



7. Conclusion 

In this paper we study the clustering of the X-ray selected AGNs in the 0.4 deg^ Chandra 
contiguous survey of the Lockman Hole Northwest region. Based on our previous study, the 
size of the CLASXS field is large enough that the cosmic variance should not be important. 
We supplement our study with the published data from the CDFN. The very similar LogN- 
LogS of the CLASXS and CDFN suggests that the cosmic variance should not be important 
when the CDFN is included in the analysis. The very deep CDFN gives a better probe of 
the correlation function at small separations. A total of 233 and 252 non-stellar sources from 
CLASXS and CDFN respectively are used in this study. We use the correlation function 
in the redshift-space as a major tool in our clustering analysis. For the whole sample, we 
have also performed an analysis using the projected correlation function. This allows us to 
quantify the effect of redshift distortion. 

We summarize our results as follows: 

1. We calculated the redshift-space correlation function for sources with 0.1 < z < 3.0 in 
both the CLASXS and CDFN fields, assuming constant clustering in comoving coordinates. 
We found a 6.7a clustering signal for pairs within s < 20 Mpc in the CLASXS field. The 

correlation function over scale of 3 Mpc< s < 200 Mpc is found to be a power-law with 
7 = 1.61q3 '^0 = 8-0ii'5 Mpc. The redshift-space correlation function for CDFN on 
scales of 1 Mpc< s < 100 Mpc is found to have similar correlation length sq = S.SSIq ^I Mpc, 
but the slope is shallower (7 = 1.3 ± 0.1). 

2. We study the projected correlation function of both CLASXS and CDFN. The best- 
fit parameters for the real-space correlation functions arc found to be tq = S.ltH Mpc, 
7 = 2.1 ± 0.5 for CLASXS field, and tq = 5.8t\i Mpc, 7 = 1.38t[!:}^ for CDFN field. Our 
result for the CDFN shows perfect agreement with the published results from Gilli et al. 
(2004). Fitting the combined data from both fields gives ro = Q-lti'^ Mpc and 7 = 1.471° ^q- 

3. Comparing the redshift- and real-space correlation function of the combined CLASXS 
and CDFN fields, we found the redshift distortion parameter (3 = 0.4 ± 0.2 at an effective 
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redshift z — 0.94. Under the assumption of ACDM cosmology, this imphes a bias parameter 
h 2.04 ± 1.02 at this redshift. 

4. We tested whether the clustering of the X-ray sources is dependent on the X-ray 
spectra in the CLASXS field. Using a hardness ratio cut at B.R — 0.7, we found no significant 
difference in clustering between hard and soft sources. This agrees with previous claims. 

5. With the large dynamic range in X-ray luminosity, we found a weak correlation 
between X-ray luminosity and clustering amplitude. Using a simple model based on ob- 
servations that hnks the AGN luminosity and halo mass, we show that the observed weak 
correlation is consistent with the model, except at low luminosities, where sources are likely 
to have lower Eddington ratio. The non-detection of a strong correlation between X-ray 
luminosity and clustering amplitude also suggests a narrow range of halo mass. 

6. We study the evolution of the clustering using the redshift-space correlation function 
in 4 redshift intervals ranging from 0.1 to 3.0. We found only a mild evolution of AGN 
clustering in both CLASXS and CDFN samples. This qualitatively agrees with the results 
based on optically selected quasars from 2dF survey. The X-ray samples, however, show a 
similar correlation amplitude as that of the 2dF sample. This is consistent with the weak 
correlation between AGN luminosity and the clustering amplitude found in this work. 

7. We estimate the evolution of bias by comparing the observed clustering amplitude 
with expectations of the linear evolution of density fluctuations. The result show that the 
bias increases rapidly with redshift (6(2; = 0.45) ~ 0.95 and h{z = 2.07) ~ 3.03 in CLASXS 
field). This agrees with the findings from 2dF. 

8. Using the bias evolution model for dark halos from Sheth, Mo & Tormen (2001), 
we estimated the characteristic mass of AGNs in each redshift interval. We found the mass 
of the dark halo changes very little with redshift. The average halo mass is found to be 

log {Mhalo/Mo) ~ 12.11. 

Our results have demonstrated that deep X-ray surveys are a very useful tool in studying 
how AGNs trace the large scale structure. Such knowledge provides an unique window to 
the understanding of AGN activity, and its relation to structure formation. The higher 
spatial density and much better completeness compared to current optical surveys allows us 
to study clustering on scales only accessible to very large optical surveys such as the 2dF 
and the SDSS. The high spatial resolution and positional accuracy of Chandra is critical for 
unambiguous optical identifications. Since our results on the evolution of AGN clustering 
could still be affected by a small number of large scale structures, as seen in Chandra Deep 
Field South, larger Chandra fields are still needed to improve the measurements. 
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Table 1. Redshift-space Correlation Function 
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Table 2. Luminosity dependance of Correlation Function 
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Table 3. Evolution of redshift-space Correlation Function 
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Table 4. Bias evolution and dark matter halo mass 
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Fig. 1. — Redshift distribution of the optically identified X-ray sources. Solid line: the 
CLASXS field; dashed line: the CDFN field. 
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Fig. 3. — The probability of source detection as a function of off-axis angle and 2-8 keV 
fluxes. Contour levels are 0.1,0.3,0.5,0.7,0.9, 0.95,0.99. Upper(lower) panels: soft (hard) 
band; Left (right) panels: 70 ks exposures and 40 ks exposures. 
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Fig. 4. — The right panels shows the random sources after detections. The pixel size is 
0.492". The left panels show the cumulated counts of the simulated sources (solid line) and 
that of the observed sources (dashed hne). Top panels: hard band; bottom panels: soft 
band. 
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Fig. 5. — The 2-8 keV flxix vs. redshift in CLASXS sample. There is no significant correlation 
between X-ray flux and redshift. 
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Fig. 6. — (a). The redshift-space correlation function for CLASXS field with 3 Mpc< 
s <200 Mpc. (b). Maximum-likelihood contour for the single power-law fit. Contour levels 
are = 2.3, 6.17, 11.8, corresponding to la, 2a and 3a confident levels for two parameter 
fits. 
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Fig. 7. — The same as Figure 6 for the CDFN except that the correlation function is calcu- 
lated for separations 1 Mpc< s <100 Mpc. 
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Fig. 8. — The projected correlation function for tlie CLASXS, CDFN fields and the com- 
bined sample, (a)-(c) show the best-fit parameters as well as the contours for the 
CLASXS-FCDFN sample, the CLASXS sample, and the CDFN sample, respectively. Con- 
tour levels are for Icr, 2(7, and Sa confident levels; (d) The projected correlation function 
for CLASXS (open circles) and CDFN (black dots) fields. Lines are the best-fit shown in 
(a)-(c). Sohd line: CLASXS+CDFN; Dotted line: CLASXS; Dashed hne: CDFN 
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Fig. 9. — Two dimensional redshift-space correlation function C(rp,7r) of the combined 
CLASXS and CDFN sample (dashed-dotted contour). Solid line shows the best- fit model. 
Both the data and model correlation functions are smoothed using a 2 x 2 boxcar to reduce 
the noise for visualization only. 
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Fig. 10. — The luminosity of X-ray sources vs. redshifts in the CLASXS (dots) and the 
CDFN (open circles) samples. 
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Fig. 11. — Luminosity dependence of clustering of AGNs. Black dots: CLASXS samples; 
Filled boxes: CDFN samples; Diamonds: 2dF sample (Groom et al. 2004). Lines are the 
models for different halo profile from Farrarese (2002). Solid line: NWF profile (k = 0.1, 
A = 1.65); Dashed line: weak lensing determined halo profile (Seljak, 2002; k, = 0.67, 
A = 1.82); Dash-dotted line: isothermal model {k = 0.027, A = 1.82) 
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Fig. 12. — The Rcdshift-space correlation function for the CLASXS field in four redshift bins. 
Left panels: The correlation functions and the power-law best-fits using maximum-likelihood 
method. Right panels: the maximum-likelihood contours for the corresponding correlation 
functions on the left. Contour levels correspond to Icr, 2a and 3(7 confident levels. 
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Fig. 13. — The Redshift-space correlation function for the CDFN field in four redshift bins, 
(layout and contour levels are the same as in Figure 12). 
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Fig. 14. — The Redshift-space correlation function for the CLASXS+CDFN field in four 
redshift bins, (layout and contour levels are the same as in Figure 12). 
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Fig. 15. — The evolution of clustering as a function of redshift for the CLASXS, CDFN and 
the two fields combined. 
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Fig. 16. — A comparison of clustering evolution in the combined Chandra fields (big dots), 
CLASXS field (big filled triangle), 2dF (diamonds), ROSAT NGP (filled box) and AERQS 
(empty box) . The solid line represent linear evolution of clustering normahzed to the AERQS. 
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Fig. 17. — The median luminosities of the 2dF quasar (COS) as a function of redshift 
(diamonds) compared to the median luminosities of CLASXS sample (triangles) and of 
CLASXS+CDFN sample (big dots). The lower panel shows the ratio of 2dF median lumi- 
nosities to the X-ray samples. 
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Fig. 18. — (a) bias evolution. The symbols have the same meaning as in Figure 16. The solid 
line is the best-fit from C05. Dash-dotted line shows the linear bias evolution model, (b). 
The mass of host halo of the X-ray sources corresponding to the bias in panel (a). 
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Fig. 19. — Evolution of Eddington ratio. Solid line: Using luminosity function from Barger 
et al. (2005). Dashed line: using luminosity function from Ueda (36) a,t z < 1.2. 



